CLAIMS 



Please amend the claims as follows: 

Claim 1 (previously presented) A computer-implemented method using a 
computer having a display and a memory for storage of information representative of 
hydrologic and hydraulic conditions for determining water flow directions based on flow 
vector values calculated in part with a D8 method applied to a watershed divided into more 
than one predefined, rectangular major cell each encompassing a minimum of nine 
predefined, rectangular minor cells including a center minor cell and eight neighbor minor 
cells wherein the watershed includes one or more types of problem areas identified as 
depressions, flat areas or peaks each of which extends across more than one major cell and 
each of which requires fixing and wherein further the number, location and elevation of all 
minor cells are known as are a maximum depth depression to be corrected, a maximum 
number of downstream minor cells to check in fixing depressions and peaks, a largest 
depression area to fix and a largest flat area to fix, comprising: 

selecting and designating one major cell for which to begin generating flow 
vector values; 

creating a buffered array in the memory containing elevation and identification data 
for all minor cells in the designated major cell and all major cells bordering the designated 
major cell; 

calculating flow vector values for each minor cell within the buffered array using 
the D8 method wherein a flow vector value of zero is assigned to any minor cell located in 
a problem area for which a flow vector value cannot be calculated using the D8 method; 

storing in the memory the calculated flow vector values in association with the 
minor cell to which they relate; 

checking the flow vector value of each minor cell within the designated major cell 
until all cells have been checked and all cells having a flow vector value of zero have been 
found; 

assigning a problem type in the memory to each minor cell within the designated 
major cell having a flow vector value of zero; 
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fixing first all minor cells in the designated major cell to which a depression 
problem type has been assigned wherein depressions are fixed independently and 
successively from the highest depression to the lowest depression and elevation changes 
are imposed on minor cells in the buffered array, as required; 

fixing second all minor cells in the designated major cell to which a flat area 
problem type has been assigned wherein each flat area is fixed independently and elevation 
changes are imposed on minor cells in the buffered array, as required; 

fixing third all minor cells in the designated major cell to which a peak problem 
type has been assigned wherein peaks are fixed independently and successively from the 
lowest peak to the highest peak and elevation changes are imposed on minor cells in the 
buffered array, as required; 

recalculating and storing in the memory flow vector values for each minor cell 
within the buffered array using the D8 method wherein a flow vector value of zero is 
assigned to any minor cell located in a problem area for which a flow vector value cannot 
be calculated using the D8 method; 

checking again the flow vector value of each minor cell within the designated 
major cell until all cells have been checked and all cells having a flow vector value of zero, 
if any, are found; 

if any minor cells having a flow vector value of zero have been found, 

assigning a problem type in the memory to each minor cell having a flow 
vector value of zero; 

fixing first all minor cells in the designated major cell to which a depression 
problem type has been assigned wherein depressions are fixed independently and 
successively from the highest depression to the lowest depression by imposing elevation 
changes on minor cells in the buffered array, as required; 

fixing second all minor cells in the designated major cell to which a flat 
area problem type has been assigned wherein each flat area is fixed independently by 
imposing elevation changes on minor cells in the buffered array, as required; 

fixing third all minor cells in the designated major cell to which a peak 
problem type has been assigned wherein peaks are fixed independently and successively 



3 



from the lowest peak to the highest peak and elevation changes are imposed on minor cells 
in the buffered array, as required; 

storing in the memory the elevation and flow vector value associated with each 
minor cell as changed, if necessary; 

selecting and designating a previously unselected major cell and repeating the 
process beginning with creating a buffered array; and 

retrieving the elevation and flow vector value for the watershed stored in the 
memory to depict on the display at least one of a basin delineation, flowpath tracing, 
drainage area calculation, the extents of hazardous waste spills and a streamline. 

Claim 2 (previously presented) The method of claim 1 wherein calculating for 
each minor cell further comprises: 

retrieving identification and elevation data for the minor cell; 

if that minor cell is located on the edge of the buffered array, assigning a flow 
vector value of zero to that minor cell; 

if a distinct flow direction from that minor cell to one of its eight neighboring 
minor cells can be ascertained, assigning a flow vector value of between one and nine to 
the minor cell wherein each number represents a different direction; and 

otherwise assigning a flow vector value of zero to the minor cell. 

Claim 3 (previously presented) The method of claim 1 wherein assigning for each 
minor cell having a zero flow vector value further comprises: 

obtaining the location and elevation of all neighbor minor cells of the minor cell; 
if the minor cell is at the edge of the buffered array, assigning the problem type as 

"edge- 
otherwise, calculating the slope from the minor cell to each of its eight surrounding 
neighbor cells in order to find a maximum slope and, 

if the elevation of the minor cell is lower than each of those neighbor minor 
cells, assigning the problem type as "depression"; 

if the maximum slope from the minor cell to one or more neighbor minor 
cells equals zero assigning the problem type as "flat"; and 
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if the maximum slope from the minor cell to more than one neighbor minor 
cell does not equal zero, assigning the problem type as "peak". 

Claim 4 (currently amended) The method of claim 1 wherein fixing first is a 
subroutine comprising: 

initializing a first variable in the memory to track a highest depression elevation 

found; 

flagging in the memory that original minor cell in the designated major cell having 
a zero flow vector value which has a highest elevation as compared with all other minor 
cells having zero flow vector values located in the designated major cell; 

retrieving from the memory location and elevation data for the flagged original 
minor cell and its neighbor minor cells; 

setting an initial depth limit in the memory to a preselected amount to limit the 
maximum depth depression to correct; 

initializing a second variable in the memory to track the number of times the 
maximum depth depression has been modified up to a preselected maximum number of 
times; 

creating a depression array in the memory mirroring the buffered array for 
identifying those minor cells which contribute to the depression; 
increasing the second variable in the memory by one; 
defining an area contributing to the depression; 
identifying potential outlets from the depression; 

determining whether any of the potential outlets from the depression is a good 

outlet; 

if no good outlet is found, 

exiting the subroutine if 

the contributing area extends beyond the buffered array; or 
the number of maximum depth modification attempts equals the 
preselected maximum number of times; 

otherwise, incrementing the depth limit in the memory by an amount equal 
to one-eighth of the initial depth limit; 
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if the depth limit is greater than the combination of the elevation of the 
flagged original minor cell and the initial depth limit, exiting the subroutine; 

otherwise, returning to creating a depression array; 
if a good outlet exists, 

choosing the good outlet having the lowest elevation; 

drawing a hypothetical straight line in the memory from the center of the 
flagged original minor cell to the good outlet cell; 

extending the hypothetical straight line in the memory to follow a 
downstream path to an elevation lower than the elevation of the flagged original minor 
cell; 

burning a cutline into the buffered array in the memory by correcting 
elevations in minor cells, as necessary; 

recalculating and storing in the memory the flow vector values for all minor 
cells in the buffered array; and 

unless all depressions are fixed or there is no more data, returning to 

flagging. 

Claim 5 (previously presented) The method of claim 4 wherein the initial depth 
limit is set at one-half of the originally specified maximum depth depression. 

Claim 6 (previously presented) The method of claim 4 wherein defining 
further comprises: 

obtaining data identifying the location and elevation of a center minor cell; 

testing each of the eight neighbor minor cells surrounding the center minor cell to 
determine if the elevation of any neighbor minor cell is greater than or equal to the 
elevation of the center minor cell; 

if so, flagging the center minor cell as a contributing minor cell and redesignating 
that neighbor minor cell as the center minor cell if the elevation of the neighbor minor cell 
is less than the initial depth limit and the neighbor minor cell is not located at the edge of 
the buffered array; 
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returning to obtaining data until all neighbor minor cells of all minor cells in the 
buffered array have been tested or until a stack space error is encountered or the center 
minor cell is determined to be located at the edge of the buffered array in which case the 
subroutine is exited. 

Claim 7 (previously presented) The method of claim 6 including the following 
after defining: 

if an error is encountered during the defining, 

exiting the subroutine if stack space has been exceeded, the area has already 
been flagged as too shallow or the preselected maximum number of times to attempt fixing 
the depression array has been reached; 

incrementing the depth limit; 

if the depth limit after incrementing is less than the elevation of the flagged 
original minor cell plus the initial depth limit, 

resetting the depression array; 

incrementing the second variable by one; and 

returning to defining; 
otherwise, exiting the subroutine. 

Claim 8 (previously presented) The method of claim 7 wherein in incrementing 
the depth limit is increased by one-eighth of the initial depth limit. 

Claim 9 (previously presented) The method of claim 4, wherein identifying 
potential outlets further comprises: 

obtaining data identifying the location and elevation of a contributing minor cell 
and of its neighbor minor cells; 

identifying which of the eight neighbor minor cells of the contributing minor cell is 
not part of the contributing area; 

testing whether the elevation of any such neighbor minor cell is less than the 
elevation of the contributing minor cell; 

if so, flagging such neighbor minor cell as an exterior cell; 
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after all neighbor minor cells have been so tested, further flagging the contributing 
minor cell as a potential outlet if at least one neighbor minor cell has been flagged as an 
exterior cell. 

Claim 10 (previously presented) The method of claim 9, wherein determining is a 
further subroutine further comprising: 

obtaining data identifying the location and elevation of a contributing minor cell 
flagged as a potential outlet and of its neighbor minor cells; 

setting a trace counter equal to zero; 

incrementing the trace counter by one; 

exiting the further subroutine if the trace counter indicates a value greater than the 
maximum number of downstream cells permitted to be checked; 

obtaining data for the location and elevation of the next downstream minor cell; 
testing whether that downstream cell is located at the edge of the buffered array; 
if so, exiting the subroutine; 

if not, further testing whether the elevation of the downstream minor cell is 
less than the elevation of the flagged original minor cell; 

if so, flagging the potential outlet as a good outlet and exiting the 
further subroutine; and 

if not, returning to incrementing the trace counter. 

Claim 1 1 (previously presented) The method of claim 1 wherein fixing third is a 
subroutine comprising: 

flagging in the memory that original minor cell in the designated major cell having 
a zero flow vector value which has a lowest elevation as compared with all other minor 
cells having zero flow vector values located in the designated major cell; 

retrieving from the memory the location and elevation data for the flagged original 
minor cell and its neighbor minor cells; 

finding and storing in the memory a maximum slope existing between the flagged 
original minor cell and one of its neighbor minor cells using the D8 method; 
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identifying each neighbor minor cell where the slope between that cell and the 
flagged original minor cell equals the maximum slope; 

marking each such neighbor minor cell in the memory as a potential path; 

if there are not multiple potential paths, exiting the subroutine; 

tracing each potential path downstream from each flagged neighbor minor cell up 
to the specified maximum number of downstream cell iterations wherein for each iteration 
the maximum slope is reset in the memory to be equal to the steepest slope between a 
downstream cell and one of its neighbor minor cells if that slope is the steepest slope found 
from that downstream cell; 

selecting a single potential path at random from all of the initial possible potential 

paths; 

designating the selected single potential path in memory as a best path if all 
potential paths have been checked and either the maximum slope for each potential path 
equals -9999 or no potential path has a slope equal to the maximum slope or there is more 
than one potential path having a slope equal to the maximum slope, otherwise, designating 
the single potential path with a slope equal to the maximum slope in the memory as a best 
path; 

lowering the elevation stored in the memory for the first neighboring minor cell of 
the flagged original minor cell along which the best path runs; 

recalculating and storing in the memory flow vector values for all minor cells in the 
buffered array; and 

unless all peaks are fixed or there is no more data, returning to flagging. 

Claim 12 (previously presented) The method of claim 1 1 wherein selecting a 
single potential path further comprises: 

identifying each potential path with a sequential number beginning at one; 

dividing the number of remaining potential paths by two to obtain a quotient; and 

selecting the potential path having the sequential number corresponding to the 
quotient. 
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Claim 13 (original) The method of claim 11 wherein the lowering step comprises 
reducing the elevation by .000001 feet. 

Claim 14 (previously presented) The method of claim 1 1 wherein lowering is 
followed by: 

verifying that lowering of the elevation does not introduce unexpected problems or 
undesirable flow direction changes; and 

if problems or changes were introduced, returning the elevation to its original level 
and exiting the subroutine. 
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